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Abstract 

The dependence of the average number of partons per clan on vir- 
tually and rapidity variables is analytically predicted in the framework 
of the Generalized Simplified Parton Shower model, based on the idea 
that clans are genuine elementary subprocesses. The obtained results 
are found to be qualitatively consistent with experimental trends. This 
study extends previous results on the behavior of the average number 
of clans in virtuality and rapidity and shows how important physical 
quantities can be calculated analytically in a model based on essentials 
of QCD allowing local violations of the energy-momentum conservation 
law, still requiring its global validity. 
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1 Introduction 



The first motivation of our search lies in the relevance of clan structure analysis 
both in theoretical and experimental study of multiparticle production. It is 
indeed remarkable that clan structure analysis of charged particle Multiplicity 
Distributions (MD's) in different reactions and at different cm. energies, in 
full phase space (fps) and in symmetric rapidity intervals, Ay, revealed an 
extraordinary series of new regularities 0, whose interpretation is still in part 
puzzling. On the theory side, it is to be pointed out how the idea of the pro- 
duction of independent intermediate objects (clans, clusters, strings...) which 
then decay following a cascading mechanism, firstly overlooked, starting from 
1986 became slowly the corner stone of any model of multiparticle production. 
The second motivation lies in the lack of a sound justification of clan parame- 
ters behavior from first principles within a parton shower model, which is the 
natural framework of the present study of multiparticle production phenom- 
ena in view of the application of Generalized Local Parton-Hadron Duality 
(GLPHD)|2j as hadronization prescription (see for instance [[|]). 

Along this line, we introduced the Generalized Simplified Parton Shower 
(GSPS) model, based on essentials of QCD and local weakening of conserva- 
tion laws, and, in this model, we calculated analytically 0] the average number 
of clans produced in a given rapidity interval Ay by an ancestor parton of 
maximum allowed virtuality W, N(Ay, W). It has been found that N(Ay, W) 
behavior predicted by the GSPS model is qualitatively consistent with experi- 
mental findings in e + e~ annihilation at hadron level. In particular, N(Ay, W) 
grows in rapidity in a way very close to linear and then bends toward a constant 
value as Ay approaches full phase space. N(Ay, W) shows also an approximate 
(5%) energy independence at ancestor energies below 100 GeV. This is actu- 
ally a very slow decrease which is reminiscent of the results obtained in Monte 
Carlo simulations for single gluon jetsj5|. This model shows also the correct 
growing of iV(fps, W) with initial virtuality W. In addition to these results, a 
new regularity for rescaled quantities is predicted: the ratio of N(Ay, W) to 
the value in full phase space iV(fps, W) is approximately energy independent 
within a much higher degree of approximation than that previously seen for 
N(Ay, W), when it is plotted as a function of rapidity interval rescaled to full 
phase space. 

These trends of the average number of clans in virtuality and rapidity 
variables are really remarkable and strongly demand to complete the research 
program announced at the end of our previous work, i.e., to calculate ana- 
lytically also the average number of partons per clan in the same variables, 
n c (Ay, W), and to check its behavior with experimental data. In order to 
do that, first we calculate analytically the average number of partons in the 
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shower as a function of rapidity interval Ay and virtuality W; then, by using 
previous findings^ on the average number of clans, N(Ay, W), we determine 
the average number of partons per clan, n c (Ay, W). 

The plan of the paper goes as follows. In Section 2 clan definition is re- 
examined in the framework of a two step process in view of the relevance of 
multiplicity channels and 1 in our analysis. In Section 3 a summary of the 
GSPS model is presented. In Section 4 the explicit analytical calculation of 
the average number of partons in the shower and per clan in the GSPS model 
is performed. Comments on the obtained results are given in Section 5. 

2 Two step processes and clan definition 

As discussed in the Introduction, the main conclusion of experimental analysis 
of clan structure parameters is that the parton production process within a 
shower is a two step process: to the initial production of independent inter- 
mediate parton sources (clans), it follows in the second step the production of 
partons inside each clan, whose average number calculation is the main subject 
of this paper. 

Following the standard procedure 0, let us assume to produce in the first 
step (independently from the second one) N objects (N = 0, 1, . . .) with prob- 
ability pn and generating function 

oo 

f(z) = £ p N z N . (1) 

N=0 

Each of the N produced objects gives origin in the second step to partons 
according to the same multiplicity distribution, q n ., (rij = 0, 1 . . ., Y&=i n i = n ), 
whose generating function is 

oo 

g(z) = £ q ni z n < . (2) 
The two step nature of the process is summarized in the equation: 

oo oo 

F{z) ee £ P n z n = £ p N [g(z)] N = f[g{z)\ , (3) 

n=0 N=0 

where P n is the probability to produce n partons and F(z) is the corresponding 
generating function. It should be pointed out that Eq. (|3|) is not limited to 
the case in which all produced clans are identical, but it requires only that all 
clans can be described by the same generating function. Suppose in fact that 
a clan's MD, q ni (£), depends on a set of parameters, denoted collectively by £: 
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if for each clan these parameters are independent of the values of other clans' 
parameters and of the number of clans generated, then it is possible to define 
an "average clan" whose MD is the MD of a single clan averaged with the 
probability distribution function (/>(£) that a clan is produced with parameters 

f QnMmn . (4) 



In this case Eq. (|3|) would still be valid with the function g(z) generating now 
the distribution q n .. This is the type of average that will be performed in the 
context of the GSPS model in Section 4.1. 

Notice that when f(z) is the generating function of a Poissonian distribu- 



tion, i.e.: 



f(z) = exp [N(z - l)j , (5) 

it follows: 



exp 



N[g(z)-1] . (6) 



These equations fully define the class of Compound Poisson Distributions 
(CPD) with all important related properties]^] and, in particular, the clan 
concept as group of particles of common ancestor ||. Notice that each clan 
contains at least one particle in order to be uniquely defined. This fact corre- 
sponds to the request that partons' MD inside a clan is shifted, i.e., q = 0. 
This remark is of particular interest for the extension of our discussion on clan 
parameters from full phase space to rapidity intervals. In fact, while in full 
phase space we are sure that each clan satisfies the condition q = 0, in rapidity 
intervals one should be careful because if all partons belonging to a clan fall 
outside the interval Ay, then q Q (Ay) ^ 0. The problem now is how to define 
clan concept in this case, i.e., how to redefine clan structure parameters with 
at least one parton inside each clan in rapidity intervals. The goal for a CPD 
is obtained by solving the equation: 

exp {iV(fps, W)[g Ay (z) - 1]} = exp {N'{Ay, W)[g' Ay (z) - 1]} , (7) 

where 

g Ay (z = 0) = q (Ay) ^ (8) 

and 

g' Ay (z = 0) ee q' (Ay) = . (9) 

Notice that the right-hand side of Eq. (0) corresponds to the standard parametriza- 
tion commonly used in fitting procedures. 
From Eq. (|7]), one has: 

N'(Ay,W) = N(ips,W)[l-q (Ay)} = -\ogP (Ay) , (10) 
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Po(Ay) being the rapidity gap probability and 

, g Ay (z) - q (Ay) 

the rescaled generating function of partons' MD inside a clan. The average 
number of partons per clan is given then by: 

,, A Ttr . n c (Ay,W) N(tpa,W) , A Ttr . 

It should be noticed that [1 — qo(Ay)] is the probability that at least one 
parton belonging to a clan falls in the rapidity interval Ay, and N'(Ay,W) 
can then be interpreted as the average number of clans in Ay. The probability 
to have N' clans in the symmetric rapidity interval Ay, pN'(Ay,W), can be 
related ||] to the corresponding probability in full phase space, pjv(fps, W), by 
means of the following equation: 

oo 

PN ,(Ay,W)= £ U(N',Ay\N,ips) PN (ip S ,W), (13) 

N=N' 

where H(N', Ay\N, fps) is the conditional probability to have N' clans in Ay 
when one has N clans in full phase space. Since in Q clans in full phase space 
are independently produced, U.(N', Ay\N, fps) turns out to be in this 
binomial distribution: 

n(iV',Ay|iV,fp S ) = (^)n(Ay,W) N '[l-n(Ay,W)} N - N ' , (14) 

where ir(Ay, W) is the conditional probability to have one clan in Ay when 
one clan is produced in full phase space. 

This result coincides with the previous one (see Eq. (|T0D); one should simply 
identify ir(Ay, W) with [1 — q (Ay)], i.e., with the probability that at least one 
parton in a clan falls in Ay. 

The discussion so far has dealt with a distribution for clans that allows 
events in which no clans are produced, i.e., empty events. While this may 
be appropriate at charged hadron level, at parton level it is not because, in 
order to have a single jet, we must have at least one parton, hence at least 
one clan (in the following we will refer to this problem as to the "contribution 
of the ancestor" because it is the ancestor parton which is always present in 
the jet). In full phase space, if additional clans are again produced according 
to a Poisson distribution, this contribution leads to the following equation for 
generating functions: 

F ips (z) = g ips (z) exp {A(fps, W)[g ips (z) - 1]} , (15) 
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with 

A(fps, W) = N{ips, W) - 1 (16) 

and 

g {ps {z = 0) = go(fps) = . (17) 
Accordingly, in rapidity intervals Ay, Eqs. (|T5l ) and (|TTD become: 

F Ay {z) = g Ay {z) exp {\{fps, W)[g Ay {z) - 1]} , (18) 

with 

g Ay {z = 0) = q {Ay) ^ . (19) 

In order to properly define clans in rapidity intervals Ay, the calculation goes 
as in the previous scheme without the contribution of the ancestor and the 
redefinition procedure gives: 

g Ay (z) exp {A(fps, W)[g Ay (z) - 1]} = exp {\'(Ay, W)[g' Ay (z) - 1]} (20) 

with g' Ay {0) = 0. 

Notice that the distribution on the r.h.s. must now allow for intervals empty 
of partons, hence it is not shifted, and also that the average number of Poisson 
distributed objects from this form results to be X'(Ay, W). It follows: 

X'(Ay, W) = A(fps, W)[l - q (Ay)} - logg (A 2/ ) = - logP (A 2/ ) , (21) 

/ fjs_i, A(fps, W)[g Ay (z) - 1] + log g Ay (z) 
9a " {Z) ~ 1 + X'(Ay, W) {22) 

and 

"c(Ay, W) = 1 n c (Ay, W) . (23) 

Equations (£T]), @ and (|23} should be compared with Eqs. QTD|) , (jT^) and 
(|T2|), where the contribution of the ancestor is not taken into account. 

It should be pointed out that the above procedure is correct only if the 
rescaled generating function of partons' MD inside a clan, g' Ay {z), is properly 
defined, i.e., if 

g\z = 1) = 1 (24) 

and 

^) « (%^_ o > . (25) 

Being g{z — 1) = 1, thanks to Eq. (^), the first condition is trivially fulfilled; 
the second condition, Eq. (p5|), can be reformulated as follows: 

A(fps, W) 1 ( d n \ogg Ay {z) \ 

q ^ y) = x^w) qn ^ y) + awv) — J^ o > o (26) 
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and it is fulfilled if 




Vn , (27) 

i.e., if log gAy(z) is a generating function up to a normalization factor. Being 
^°E,(gAy{z)/q (Ay)) the generating function of combinants associated in Ay 
with partons' MD inside a clanj7|, gA y ( z ), the validity of Eq. (|27|) for all n 
implies that combinants are positive definite, i.e., that the corresponding MD 
inside a clan is also a CPD 0. 

Notice that now the binomial convolution of Eq. (|T3| ) , applied in [|J , gives 
a result different from Eq. (pl|): 



N"(Ay, W) = [l- q (Ay)}N(fps, W) ± X'(Ay, W) ; (28) 

it is N"(Ay, W) which should properly be identified as the average number 
of clans that produce at least one parton in the interval Ay. Thus in the 
present case in which at least one clan is always produced one cannot simply 
identify the parameter which appears in a CPD description of the MD in a 
rapidity interval with the average number of clans actually contributing to 
that interval: the former, given in Eq. (pl|) diverges when the interval Ay 
goes to fps; the latter, given in Eq. (p8|), tends to iV(fps, W) in the same 
limit. This point does not arise in the experimental measure of charged hadron 
MD for the simple reason that events with no charged particles in full phase 
space are allowed (although they are not measured). The detailed analysis 
we are carrying out in this paper has brought this problem to the surface: 



earlier work on the MD at parton level |10[ has by-passed this difficulty by 
neglecting multiplicities and 1 in full phase space: this was indeed justified 
as the interest was mainly in small rapidity intervals and, more important, 
the concept of clans was still a statistical one. But in making the step and 
treating clans as genuine subprocesses, this resolution is no longer appropriate; 
the consequences of this fact will be discussed in Section 5 together with our 
results. 

Having clarified this important point on clan definition in a two step pro- 
cess, we can now proceed to calculate in the framework of the GSPS model 
the average number of particles inside a clan as a function of rapidity and 
virtuality variables, i.e.: 



n c (Ay,W) = d9A ^ 



dz 



(29) 



It should be clear that clan parameters used in experimental analysis differ 
from those calculated in the GSPS model via Eq. (p9|). In fact, in experimental 
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analysis, one determines g' Ay (z) (Eq. (p2[)) and, accordingly, 



<(Ay,W) 



dz 



z=l 



(30) 



The primed quantities of Eq. ([!(]) are of course linked to the unprimed ones 
of Eq. © by Eqs. (0)- (||). 

It is to be pointed out that the average number of particles in a given 
rapidity interval Ay, n(Ay, W), can be determined both before and after the 
redefinition process. One has: 



n(Ay, W) = A(fps, W)n c (Ay, W) 



(31) 



and 



n(Ay,W) = X'(Ay,W)n' c (Ay,W) 



(32) 



respectively; notice that, as expected, n(Ay, W) is independent of the redefi- 
nition procedure. 



3 The GSPS model: a summary 

The description of a parton shower initiated by an ancestor parton of virtuality 
Q splitting into two partons of virtualities Qo and Q\ in terms of virtuality and 
rapidity variables is usually given within LLA QCD evolution equations[ll] by: 

dP Q^QoQi _ a 8 (t) T p , w /oo^ 
-j t ~2~ J ^Q-^QoQi\ z a) az o i K 66 ) 

^-min 

where 

12vr , Q 2 

as{t)= (llN c -2N f )t > t = l ° g A^ ■ (34) 

N c and Nf are the number of colors and flavors respectively and A is the QCD 
scale. 

Pq->QoQi ( z ) i s the splitting kernel and z is the energy fraction carried away 
by the parton of virtuality Qo. The virtuality at which the parent parton 
splits as well as the energy fraction carried away by the daughter partons is 
determined at each splitting. The main problem in the standard approach 
concerns how to take care of energy-momentum conservation law (zq + z\ — 
1). In the GSPS model, a different approach is followed. The branching 



process is described here in fact in terms of finite splitting functions fl2|1 and 
this procedure allows to fix kinematical limits on the energy fraction z$ by 
two-body kinematics 1 . 
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For the virtuality dependence of the splitting function, a factorized form 
inspired to QCD and suitably normalized by a Sudakov form factor term has 
been chosen: 



PA(Qo\Q)dQ = P$(Q)C A (Q)dQ = Ca(Q)^ [c^Q^) dQo = d ( bj^r) ' 

° g (35) 

where A at this level of investigation is the only free parameter of the model 
and controls the length of the cascade. 

We include next clans as independently produced intermediate parton sources 
and modify accordingly the kinematics of the process by allowing local vi- 
olations of the energy-momentum conservation law, still requiring its global 
validity, i.e., offspring partons of virtualities Qi can fluctuate according to: 

Qo + Qi^Q , 1 GeV <Qi<Q [i=0,l] . (36) 

Under this assumption, the joint probability density V(QoQi\Q) becomes fac- 
torized: 

V(Qo,Qi\Q)dQ dQ 1 =p A (Q \Q)dQ p A (Q 1 \Q)dQ 1 . (37) 

Each parton emission is therefore independent in virtuality. Of course, con- 
straints on the energy fraction carried away by daughter partons are also no 
longer valid, i.e., zq + Z\ ^ 1, and kinematical bounds in rapidity are loosened: 

\yi-y\<\og^ [i=o,i]. (38) 

This new condition modifies the splitting kernel in z which now is decoupled: 

P(z , zi)dzodz\ oc — - — - . (39) 

Zq Z\ 

By combining Eqs. (|37D and (|3~9"D, one sees that the total bi-dimensional split- 
ting function is also decoupled and corresponds to two independent emissions 
of a single parton, i.e., in order to describe the production process, it is enough 
to follow just one branch of the shower evolution. Eq. (|3^) can be rewritten 
in terms of rapidity variables for each branch: 

dzi 

Y(\yi - y i+1 \,Qi, Q i+1 )dyi oc — = dy { , (40) 

Zi 

where Y(\yi — y i+ i\,Qi,Q i+ i) is the probability that a parent parton with 
rapidity yi + i generates in a single step a parton of rapidity yi, being their 
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virtualities Qi+i and Qi respectively. When properly normalized, Eq. (|40| ) 
turns out to be: 



dy% 

Y{\yi - y i+1 \,Qi,Qi +1 )dyi = ' r 9{\og(Q i+1 /Qi) - \yi - y i+1 \) . (41) 

21 °gw) 

By identifying now in each splitting one of the produced partons as the ancestor 
parton of a clan, one obtains that clans' production process is Markoffian and 
therefore clans are independently produced; notice that this a property of the 
model and not an assumption, as discussed in M. The explicit calculation of 
clans' MD in a parton shower originated by an ancestor parton of virtuality 
W gives in fact a shifted poissonian distribution with average number 

N(fps,W) = l + A\ogU^j . (42) 

Results of this scheme of parton shower evolution have been summarized in 
the Introduction. 

We concluded that the observed qualitative agreement of the predictions 
of the model with experimental behavior for the average number of clans in 
rapidity intervals supports our idea that clans can be considered independent 
intermediate parton sources, i.e., that local fluctuations in the virtualities of 
produced partons can occur along the shower. This fact is quite acceptable 
since energy-momentum conservation laws are expected to weaken in limited 
rapidity intervals whereas their role is surely fundamental in full phase space. 

Finally notice also that the approximation used in the GSPS model closely 
resembles that used in the dipole model [O, where two subsequent dipole emis- 
sions are needed in order to describe a full branching. Furthermore, the kine- 
matical structure of the model is similar. The main difference among the two 
models lies in the description of the shower evolution: in fact, it is based on 
finite splitting functions in the GSPS model, whereas it follows the standard 
approach with elementary splitting kernels in the dipole model. 



4 The average number of partons per clan in 
the GSPS model 



4.1 The structure of the calculation 

In order to study the average number of clans in the rapidity interval Ay, 
N(Ay, W), in (4| we limited our discussion to the first step of parton shower 
evolution in the GSPS model. It is clear that if one wants to calculate the 
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average number of partons per clan in the same interval, n c (Ay, W), one has 
to analyze the second step of parton shower evolution, i.e., to study the pro- 
duction of partons inside clans. In order to do that, inspired by the criterium 
of simplicity and previous findings flOl, |14|| , we decided to maintain inside a 



clan the structure of the model seen in the first step. The only difference (see 
Eq. (|35|) ) lies in the introduction of a new parameter, a, controlling the length 
of the cascade inside a clan, in the expression of the probability that a parton 
of virtuality Q emits a daughter parton in the virtuality range [Q , Q + dQ ], 
i-e-, 

'logQoV 



Pa(Qo\Q)dQ = d 



logQ 



(43) 



The model becomes of course a two-parameter model (A and a controlling the 
length of the cascade in step 1 and 2 respectively). 

Let us now consider a clan splitting at virtuality Q and study the evolution 
equation for the generating function of partons' MD in full phase space inside 
a clan, g{ ps {z,Q), which, according to the pure-birth structure of the model, 
turns out to be: 



dg{ ps (z,Q) 
dQ 



Pa(Q\Q) 9ip S {z, Q) - #fps(z, Q) 



where p a (Q\Q), from Eq. fl43|), is given by 

Pa(Q\Q) 



QlogQ 



(44) 



(45) 



Notice that, as expected, Eq. (fP|) corresponds to the evolution equation for the 
generating function of gluons' MD in a gluon jet in Leading Log Approximation 
of QCD with fixed cutoff infrared regularizationjTSJ, as well as to the limit 



of the Simplified Parton Shower model when one weakens conservation laws 
according to Eq. (I 



ID | . Of course Eq. ([0]) concerns now partons' MD inside 



a single clan. 

The solution of Eq. 

9f ps (z, Q) = 



is 



z 



l + (n(Q)-l)(l-*) 

z 



where 



n(Q) 



3 A a (Q) 



\a(Q) 



Q 



Q>2 GeV 
Q < 2 GeV , 

p a (Q'\Q')dQ' ■ 



(46) 



(47) 



The solutions correspond to a shifted geometric distribution (Q > 2 GeV) 
and to a clan with only one parton (Q < 2 GeV). The bound 2 GeV is a 
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consequence of the fact that in the GSPS model the virtuality cutoff is fixed 
at 1 GeV (a parton with virtuality Q < 2 GeV cannot split any further by 
assumption). Notice that this finding agrees with the clan model discussed in 



14] , where logarithmic MD for partons inside average clans is interpreted as 
the result of an average over geometrically distributed single clans of different 
multiplicity. 

Equation fl4"E| ) answers our first question and gives the generating function 
of partons' MD inside a single clan in full phase space. In order to match 
our theoretical calculations with experiments, one should introduce the rapid- 
ity dependence in the model and average single clan properties over clans of 
different rapidity and virtuality. According to our idea that single clans are 
independent intermediate parton (gluon) sources acting both in virtuality and 
rapidity spaces, and not only statistical objects as in the original clan model, 
we decided first to extend the model defined in full phase space to rapidity 
intervals and then to average these results on single clans' MD's. 

Following the discussion contained in Section 2, we calculate the generat- 
ing function of partons' MD inside a clan in rapidity intervals Ay through a 
binomial convolution on the corresponding generating function in full phase 
space: 

f n \ l + (z- l)7r B (Ay, Q, y) 

gAy[z,Q,y) = — , \ f -fn\ — TYTi V Q > 2 GeV ; 

1 + 7r a (Ay, Q, y){n{Q) - 1)(1 - z) 

= l + (z-l)ir a (Ay,Q,y) Q < 2 GeV , (48) 

where n a (Ay,Q,y) is the probability that a clan of initial virtuality Q and 
rapidity y produces a daughter parton inside the interval Ay. Notice that 

7r a (Ay,Q,y) = [ p(y f \Qy)dy f . (49) 

J Ay 

PilJflQy) i s the conditional probability that a final parton produced by a clan 
of given Q and y has rapidity yj. 

We have now extended formally the generating function of partons' MD 
inside a clan from full phase space to rapidity intervals for a single clan of 
known initial virtuality Q and rapidity y. It is clear that in this model a single 
event contains many clans of different initial parameters, while at experimental 
level one measures properties of average clans, in the sense indicated in Eq. (H). 
Accordingly, we proceed now to average Eq. ( fj8|) over Q and y. It follows that: 

rW roo 

9Ay(z,W)= dQ dyg Ay (z,Q,y)a(Qy\W), (50) 

Jl J — oo 

where gAy{z, W) is the generating function of partons' MD in the interval Ay 
for an average clan generated in a shower of virtuality W; a{Qy\W) is the 
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probability that a parton of maximum allowed virtuality W produces a clan 
of virtuality Q and rapidity y and gA y (z, Qy) is given by Eq. ( ]48|) . 

By deriving then Eq. (pjQP with respect to 2 and choosing -2=1, one gets 
the average number of partons in an average clan generated in a shower of 
virtuality W in the rapidity interval Ay, i.e.: 

rW roo 

n c (Ay, W)= dQ dy n c (Ay, Q, y)a(Q, y\W) , (51) 
J 1 J — 00 

where n c (Ay, Q,y) is the parton average multiplicity in a rapidity interval Ay 
for a single clan of given Q and ?/, and is related to n c (fps, Q, y) as follows: 

n c (Ay, Q, y) = n a (Ay, Q, y)n c (fps, Q, y) = 7r a (Ay, Q, y)e x( - Q) . (52) 

Equation (^1|) formally solves our main problem; the final analytical expres- 
sion for n c (Ay, W) results after having determined explicitly the two factors 
n c (Ay, Q, y) and cr(Q, y\W). Notice that the integration domain in y is deter- 
mined by kinematical constraints in a(Q,y\W), which are explicitly discussed 
in the following Subsection. 

In Subsections 4.2 and 4.3, we determine a(Q,y\W) and n c (Ay, Q, y) re- 
spectively, whereas in Subsection 4.4 we perform the final integration over Q 
and y according to Eq. (|5"T[). 



4.2 The probability to produce a clan of given virtuality 
Q and rapidity y, a(Q,y\W) 

The weight function a(Q,y\W) we are looking for coincides with the bi- 
dimensional clan density in virtuality Q and rapidity y for a shower origi- 
nated by a parton of maximum allowed virtuality W, Px(Q,y\W), up to the 
normalization factor iV(fps, W): 

*{Q,V\W)= N{fpStW) ■ (53) 

Notice that Ps(Q, y\W) by definition depends on the first step of the production 
process and therefore the expression of a{Q,y\W) is expected to contain the 
parameter A only. 

Following the procedure used in ||], where we calculated the clan density 
in rapidity, we define here Pz(Q, y\W) as: 

00 

MQ,y\w)= Epn(Q,v\w), (54) 

N=l 
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where Pn{Q, y\W) is the probability that a parton of maximum allowed virtu- 
ally W generates a clan of virtuality Q and rapidity y after N steps. 
The explicit formula for Pn(Q,u\W) in the GSPS model is given by: 



/max{2,Q} dQ N _ 2 PA (Qn-1 I Qn^) j dy N _ 2 Y( \y N ~2 ~ yN-1 \ , Qn~2, Qn-1 ) 



Following the GSPS model prescriptions, Y(\yi — y i+1 \,Qi,Q i+1 ) can be ap- 
proximated by a gaussian function, i.e.: 



Y(\Vi - Vi+i\,Qi,Qi+i) = 




This approximation allows to perform the integration over all rapidity variables 
with the exception of t/jv-i- The integrations on the corresponding virtuality 
variables can be done by using the symmetry properties of the integrand. One 
gets: 



pn(Qv\w) 



rW 

Jmax{2,Q} 




(55) 



Jmax{2,Q} d Q iP a{Q i\Q 2) / d yi Y(\ yi - y 2 \ , Qx, Q 2 ) 

J—oo 

Pa(Q\Qi) Y{\y- yi \,Q,Q x ). 



p N (Qy\W) = 



JQ J-oc 



rW roo 1 

/ dQ N ^p A {Q N ^\W) / dy N ^S(y N ^ - tanh-yi - (Qn^/W) 2 ) 




for Q > 2 GeV 



(57) 



and 



p N (Qy\W) = 




(58) 
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Notice that the last integration over i/n-i needs not to be performed by using 
the gaussian approximation; in fact, the integration can be done with Y(\yi — 
y i+ i\,Qi,Q i+ i) given in its original form (Eq. ((41])). This fact restores global 
conservation laws in rapidity space. 

By using Eq. @ for p A (Qo\Q) and Eq. (g7|) for X A (Q), Eq. @ leads to 
the following expression for pz(Q, y\W)\ 

MQv\w) = P$(Q)C A (W) / dQ N - X i4(Q N - X )C A (Q) 

dy N -i5(y N „i - tanh" 1 ^! - (Q N _ x /W) 2 )Y{\y - y N „ x \,Q, Q N - X ) 



+ p A (Q\W)5(y - tanh-yi - (Q/W) 2 ) 

for Q > 2 GeV (59) 



and 



MQv\W) = p£(Q)C A (W)J dQ N - X p£(Q N - X )C A (2) 

/oo , 
dy N - X 5(y N - X - tanh-Vl - (Q N „ x /W) 2 )Y(\y - y N -i\,Q, Q 
-oo 



+ p A (Q\W)5(y - tanh-yi - (Q/W) 2 ) 

for Q < 2 GeV . (60) 



By integrating over yjsr-i, one has: 



MQV\W) = p A (Q\W)C A (m a x(2,Q)) / dQ N „ x 



max(2,Q) 2 log Qn-i/Q 



9(logQ N - X /Q -\y- tanh-yi - (Q N - X /W) 2 \) 
+ p A (Q\ W)8(y - tanh-yi - (Q N - X /W) 2 ) . (61) 

Equation (|BT]) in the approximation 



tanh-yi - (Q N - X /W) 2 ~ log(2iy) - logQ^-i , (62) 

gives: 

/•log IT 

Pe(Qi/|W) = p A (Q\W)C A (Q) d(\ogQ')Q> 

Jl/2{log2W-y+logQ 2logQ'/Q 



, p£(Q') 



+ p A (Q\W)6(y - tanh-yi - (Q'/VT) 2 ) 
for Q > 2 GeV (63) 
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and 

MQv\W) = PA (Q\W)C A (2) 



I , , ^ of" Vv /n ^K^/Q ~ \V ~ tanh-yi - (Q'/Wy\) 

7max(2,Q) 2i0gQ'/Q V 



+ Pa(Q|W0% - tanh-yi - (Q'/Wy) 

for Q < 2 GeV . (64) 

The region where pz{Qy\W) is different from zero, according to the just men- 
tioned approximation, is given by 

W 2W 
- log — < y < log — . (65) 

These results allow to calculate a(Q,y\W) from Eq. (^). 

4.3 The average number of partons in Ay inside a single 
clan, n c (Ay,Q,y) 



In view of Eqs. (|49| ) and (p2]), one has now to calculate the conditional proba- 
bility that a final parton produced by a clan of given virtuality Q and rapidity 
y has rapidity p(yf\Qy). This calculation can be performed in terms of the 
probability that a clan of given virtuality Q and rapidity y produces a final 
parton of rapidity yj after iV steps, T^{yf\Qy). It follows: 

oo 

p(yf\Qy) = Yl r N{yAQy) • (66) 

N=l 

Notice that r N (yj\Qy) is obtained by integrating r N {Q fyf\Qy) over the virtu- 
ality of the final parton, Qf in the interval 1 GeV -=-2 GeV: 

rn{yf\Qy) = J dQ f r N (Q f y f \Qy) . (67) 

t n{Q 'fVflQy) is i n f &c t t ne probability that a final parton of virtuality Qf and 
rapidity yj is generated by a clan of virtuality Q and rapidity y after N steps. 
The problem is now how to calculate r^(Q fyf\Qy)- This calculation can be 
performed in the GSPS model; in fact, being 



DC 



rN(QfVf\Qy) = i dQ N ^ 1 p a (Q N ^ 1 \Q) / dy N ^xY{\ 

yN—i ~ y\i Qn—ij Q) 

J —oo 

(68 

l POO 

dQiPa(Qi\Q2) / dy 1 Y{\y 1 -y 2 \iQiiQ2) 

J —oo 

p a (Qf\Qi)Y(\y f - yi\,Qj,Qi) , 



16 



by using the same technique discussed in the previous Subsection (Gaussian 
approximation of Y(\yi — Vi+i\, Qi, Qi+i)), one finds: 

r N (Q fVf \Qy) = p a (Q f \Q)Y(\y f - y\,Q f , Q) ^ . (69) 

The structure of this calculation is very similar to that performed in 0] (Eq. (48)). 
The main difference lies in the fact that here we are studying the production 
mechanism for partons inside a clan and not for clans themselves, i.e., the 
result must depend on the new parameter a only. 

Finally, one should notice that in this scheme the ancestor parton after the 
splitting continues on its way whereas the daughter parton gives origin to the 
clan. Therefore partons' MD inside the clan is shifted by one unit. 

Our calculation can be simplified by assuming Qf — 1 GeV in Y(\yf — 
y\,Qf,Q), an assumption which is verified for clans of high virtuality and is 
quite acceptable in general. Accordingly, one has: 

p(Vf\Qv) = e XaiQ) [ dQ fPa (Qf\Q)Y{\y f - y\,Q f , Q) (70) 
which, by using Eqs. (BSf) and Q4"T|), becomes 



rivAQv) = ^q 9 ^Q -\vf-v\)- ( 71 ) 

Equation ([FTP is valid for Q > 2 GeV. For Q < 2 GeV, one produces clans 
with only one parton, i.e.: 

P Q <2(y f \Qy) = 6{y-v f ). (72) 

Having found PQ<2(yf\Qy), we can now calculate from Eqs. ([52]) and (|49|) the 
average number of partons in a single clan: 



n c (Ay,Q,y) = e Aa(<3) max 



min(sup(Ay),?/ + logQ) - max(inf(Ay), y - logQ) 



and 



21ogQ 

for Q > 2 GeV (73) 



n c {Ay,Q,y) = e A " (Q) y E Ay 
= y^Ay 
for Q < 2 GeV (74) 



respectively. 
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4.4 The average number of partons in Ay in an average 
clan 



The main content of this Section concerns the explicit analytical expression of 
n c (Ay, W). This result follows by solving the double integral over Q and y in 
Eq. n c (Ay,Q,y) is given here by Eqs. ( f73[) and flT4|), whereas a(Q,y\W) 
is obtained by inserting Eqs. (|6~3| ) and ([64]) into Eq. (|53l). 

Attention should be paid to the integration domain of y and Q, which 
is determined by the conditions 1 GeV < Q < W and —\ogW/2Q < y < 
\og2W/Q (see Eq. (§§)). Notice that according to Eqs. fl7|) and Q7§ the 
integrand has a different structure in the intervals 1 GeV < Q < 2 GeV and 2 
GeV < Q < W GeV. Accordingly, we perform separately the two integrations: 



n c (Ay, W) = n c (Q < 2) + n c (Q > 2) 
In the interval 1 GeV < Q < 2 GeV, one has: 

Aog2W/Q 



n c (Q < 2) 



dQ 

3£ W 



logW/2Q 



dy6(y c -\y\) PA (Q\W)C A (2) 



(75) 



(76) 



dlogQ'Q' P ^ Q '] + 

max[lo g 2,(lo g 2Ty-i/+logQ)/2] 2 logQ'/Q 



+ 



dQ 



\og2W/Q 



dy6(y c - \y\)p A {Q\W)8{y - tanlT 1 ^ - (Q/WO 2 ) . 



-logiy/2Q 

In the approximation of Eq. (|62"D, Eq. (|76|) can be rewritten as follows: 



n c {Q < 2) 



dQp A {Q\W)C A {2)Y, A {Q,y,W) + 

A 



(77) 



+ 0(y c -logW) 



log 2 
logW 



log2W-y c 
logVT 



where 



s A (Q,y,iy)= / 



lo R W 



dhgQQ — — — 

(lo g 2H/-y+logQ)/2 2log Q'/Q 



-(logg)^- 1 hogx + Y: 

1 I Jl=l 

In particular, for A=2, one gets: 

'hgW 



A - l\ x r 



n 



n 



logW/logQ-l 
(logW-j/+LogQ)/21ogQ-l 



(7f 



Z\Q,y,W) = logQ 



log 



logQ 



- 1 - log 



'log2W - y + \ogQ 
, 2 logQ 



- 1 



+ ! + ilogW72Q 



(79) 
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Coming to the interval 2 GeV < Q < W, one can split the integration in 
two different terms: 



n c (Q > 2) = n c (Q > 2; a) + n c {Q > 2; /?) , 



*0) 



where 



r w r 
n c {Q > 2; a) = J dQ J 



r \og2W/Q e Aa(Q) 

dy 



max 



-logW/2Q iV(fpS, W) 

min(sup(Ay), y + log Q) — max(inf(Ay), y — log Q) 



2\ogQ 



p A (Q\W)S(y -t^nh 



-i 



(Q/wy 



(81) 



and 



n c (Q>2;/?) 



ir 



dQ 



log ( 

dy 



Xa(Q) 



max 



logW/2Q iV(fps, W) 

min(sup(A?/), y + log Q) — max(inf (Ay), y — log Q) 



21ogQ 



PA(Q|^)c j4 (g)E A (g, 2/ ,w) 



(82) 



So far the rapidity interval has been generic, but in order to continue our 
explicit calculations in agreement with standard practice 0, we specialize in 
the following to central symmetric intervals of rapidity of half-width y c : Ay = 

[-VcVc]- 

By using Eq. (|62f) and the properties of the 5-function, one has: 



n c (Q > 2; a) 



dQ: 



C a (2)p£(Q)C A (W) 



21o g g C a (Q)N(ips,W) 
max [0, mm(y c , log 2W) — max(— y c , log 2W — 2 log Q)} . 



3) 



Notice that the expression of the maximum in Eq. ( p3[ ) assumes the following 
different values in different ranges of logQ, i.e., for logQ < |(log2W / — y c ), 
y c - log 21^ + 2 log Q for |(log 2W - y c ) < log Q < |(log 2W + y c ) and 2y c for 
| (log 2W + y c ) < log Q < log 2W. It follows: 



n c (Q >2;a) = 
'\og2W + y c 



1 



(log2) a (logW) A N(W) 

, A+a—l 



(y c -log2W) 



A 



2(A + a- 1) 



log2W-y c 



A+a-l 



+ 



+ 



A 
A + a 



'\og2W + y c 



A+a 



'\og2W-y c 



A+a 



+ 



54) 
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+ Vc 



A 



A + a-1 
for y c < log Wf 2 and 

n c (Q > 2; a) 



(log 



A+a-1 



'log2t^ + ?/ c 



A+a-l' 



1 



(\og2) a (\ogW) A N(W) 



{y c -\og2W) 



A 



2(4 + a - 1) 



+ 



(iogiy) A+a " 1 -(io g 2) A+a - 1 

A 



A + a 



(logW0 A+a -(log2) 



A+a 



for y c > logW/2. 

Now, we calculate the second term, n c {Q > 2;/3). It is convenient to 
introduce a new set of variables, i.e.: 

£ = logQ + y (86) 
V = logQ-y. (87) 

In these variables the initial integration domain turns out to be: 

W 

V < log( y ) , £ < log(2W) , I/ + f > 2 log 2 , (88) 
as indicated in Figure ([[]) (solid line). 

Notice also that the expression of the maximum contained in Eq. ( ff3|) assumes 
different values in different sectors of the integration domain. 

For y c in the interval (log 2 +- logW/2) one can identify four sectors (see 
Figure (|1|), dotted lines): 
I) 



7/ c <£<log2 y c < 7/ < logW/2 ; 
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II) 



-y c + 21og2 < f < t/ c t/ c < 7/ < logW/2 (90) 

max(-7/ c , - log W/8) < £ < -y c + 2 log 2 - £ + 2 log 2 < 77 < log W/2 ; 



///) 



7/ c <£<log2iy -7/ c + 21og2 < 77 < y c (91) 

-77 + 2 log 2 < £ < log 2W max(-7/ c , - log W/2) < 77 < -y c + 2 log 2 ; 



/V) 



7/ c + 21og2 < £ < 7/ c 



£ + 21og2 <r]<y c 



(92) 
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where n a (Ay,Q,y) is: 



I) 2y c (93) 

II) i + Vc (94) 
///) V + Vc (95) 
IV) i + V (96) 



respectively. 

For y c values in the extreme allowed interval, y c > \ogW/2, sectors /) and 
II) disappear and the other two sectors III and IV are modified as follows: 

HI) 

y c <£<\og2W - e + 21og2 < v < \ogW/2 ; (97) 

IV) 

-logW/8 < £ <y c -£ + 21og2 < rj < \ogW/2 . (98) 

Finally, for y c < log 2, sector IV) only disappears and the others sectors 
become: 

I) 

y c < £ <\og2W -y c + 21og2 < rj < \ogW/2 (99) 
-r] + 2 \og2 <\og2W y c < rj < -y c + 2 log2 ; 

II) 

-y c <i<y c -e + 21og2<r ? <log^/2; (100) 

HI) 

-r 1 + 2\og2<i<\og2W - y c < rj < y c . (101) 

The final analytical result for n c (Q > 2; f3) can then be obtained. In order to 
have an impression of the structure of this formula, we show it explicitly in 
the domain log 2 < y c < log W: 

UciQ 1 + Alog(^) \' Vc + + 12(1^2) + 

29 (log 2) y c 2(log2) 2 y c 7y c 3 | 

12 (log W) 2 3(logW^) 2 12 (log Wf (log 2) 
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4(logW02/ c log((logWQ) 2y c log((logWQ) 



3 (log 2)' 



(log 2) 



+ 



+ 



2 (-2 (logWQ 3 ?/ c + 3 (log HQ 2 (log 2) y c - (log2) 3 y c ) log((log HQ - (log 2)) 

3 (log W/) 2 (log 2) 2 



(-{\ogW)+y c f log((logW0-y c 
6(logiy) 2 (log2) 2 



-{\ogWf -A{\ogW) (log 2) 



6 (log 2) 2 + 2 (log W ) y c + 4 (log 2) y c - y c 



+ 



((log HQ - (log 2) - y c f log((log WQ - (log 2) - y c ) 



+ 



12 (log Wy (log 2) 
(-(log - (log 2) + y c f MOogMQ + (log 2) - y c 
12 (log W/) 2 (log 2) 2 



2 (log WO 2 + 



+ 15 (log W) (log 2) - 5 (log 2) 2 - 4 (log W) y c - 3 (log 2) y c + 2 y c 



+ 



+ 



{{\ogW)+y c f log((logWQ + y c ) 
6 (log Wf (log 2) 2 



(logiy) 2 + 4(logIV) (log 2) + 



+ 6 (log 2) 2 + 2 (log W) y c + 4 (log 2) y c + y c 



+ 



HlogjQ + (log 2) - y c y log((log - (log 2) + y c ) 



+ 



12 (log Wf (log 2) 
((log WO + (log 2) + y c ) 2 log((log + (log 2) + y c 



12 (log WO 2 (log 2) 2 



-2 (}ogW) 2 + 



15 (log W) (log 2) + 5 (log 2) 2 - 4 (log W) y c - 3 (log 2) y c - 2 y c 



+ 



(102) 



Having now integrated, as proposed in Eq. fl5lD , the product of the average 
parton multiplicity in a rapidity interval Ay for a single clan of given virtuality 
Q and rapidity y, times the probability that a parton of maximum allowed 
virtuality W produces a clan of virtuality Q and rapidity y, by dividing the 
integration domain in several regions, one obtains the wanted final analytical 
expression for the average number of partons inside an average clan in the 
interval Ay by assembling the various pieces as follows: 



n c {Ay,W) = n c (Q < 2) + n c {Q > 2) 

= n c (Q < 2) + n c (Q > 2; a) + n c (Q > 2; (3) 



(103) 



The average number of partons in the shower, n(Ay, W), can then be calcu- 
lated from Eq. (|31"|). 
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5 Results: comments and discussion 



The motivation of this work - as pointed out in the Introduction - has been 
to justify from first principle clan structure analysis in multiparticle produc- 
tion and to provide a sound theoretical basis to the regularities discovered by 
its application to experimental data. In order to do that, one would need 
QCD calculations of a bi- dimensional shower evolution in a region where QCD 
predictions are not available and calculations are very hard to be done. One 
should therefore fully rely on models. Our program was to study analytically 
clan structure in rapidity and virtuality variables in a parton shower model 
based on essentials of QCD in a correct kinematical framework, the GSPS 
model. Results of analytical calculations of the average number of clans be- 
havior in virtuality and rapidity variables were reported in M and summarized 
in the introduction. In order to complete our initial program we propose in 
this paper to discuss the properties of the second parameter of clan struc- 
ture analysis, i.e., of the average number of partons per clan, in rapidity and 
virtuality variables as resulting from the analytical calculations presented in 
previous Sections. 

Notice that these calculations have been performed by following very mild 
mathematical assumptions and by choosing for the parameters A and a of 
the GSPS model the values 2 and 1 respectively. It should be remembered 
that A and a control the "length" of the parton shower and of clans in the 
shower. The suggested choices avoid nasty inessential calculations and make 
possible the analytical solution of the model without hurting its logic. It 
should be pointed out that MC calculations based on the same architecture of 
the GSPS model allow to predict the behavior of n c (Ay, W) as well as that of 
the average number of clans in the same variables, N(Ay, W), for all values of 
the parameters a and A. In addition MC calculations here can be considered 
as exact , i.e., the result is obtained without any of the above mentioned mild 
mathematical approximations used in the analytical calculation in order to get 
the explicit dependence of N(Ay, W) and n c (Ay, W) on Ay and W . This fact 
should be considered important for estimating the goodness of the assumptions 
performed in the analytical calculations. 

Knowing from previous work[|J the rapidity and virtuality dependence of 
the average number of clan, the analytical calculation of the average number 
of partons per clan in the same variables can be done in terms of the aver- 
age number of partons in the full shower. This is indeed a quantity which - 
as previously discussed - does not depend on the redefinition procedure and 
therefore can safely be used to our purpose. 

The result of the analytical calculation of the average number of partons 
in the shower, n(Ay,W), as a function of the rapidity interval Ay and of 
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the maximum allowed virtuality W with A=2 and a=l, Eq. (|29|), is shown in 
Figure (D for W= 50 GeV (solid line), W=100 GeV (dashed line) and W= 
500 GeV (dotted line). The average parton multiplicity grows almost linearly 
with rapidity for relatively small Ay intervals and then it is slowly bending for 
Ay intervals approaching to fps, where it reaches its maximum (n(fps) ~ 10 
for W = 50 GeV, n(fps) ~ 12 for W = 100 GeV and n(fps) ~ 17 for W = 500 
GeV). 

Monte Carlo calculations of n(y c , W) for the same W values and for the 
same choice of the parameters A and a are given in Figure (|3|). The general 
trend is the same as that in Figure (|2]). The agreement between the two calcu- 
lations is very good in fps. For smaller Ay intervals Monte Carlo calculations 
are sistematically higher into few percent and allow to estimate the amount of 
the approximations done in the analytical calculation. 

It is interesting to remark that the normalized average number of partons 
in the shower, n(Ay,W)/n(fps,W) scales in virtuality as a function of the 
rescaled rapidity interval, Ay/fps, see Figure (|). This scaling in W is found 
to depend on the parameter a, as different values of a give different scaling 
curves, differently from the scaling found in || for N(Ay, W)/N(ips, W) which 
is independent of the mechanism inside clans. 

Figure (|5|a) shows the behavior of the average number of clans, N'(Ay, W), 
in showers of maximum allowed virtuality W equal to 50 GeV (solid line), 100 
GeV (dashed line) and 500 GeV (dotted line) as a function of rapidity interval 
Ay (Eq. (^)). The analytical calculation has been done with A=2 . Being 
the first step only of the production process involved in the calculation, the 
result does not depend of course on the parameter a. The Figure content is 
the same discussed in Figure 8 of ||. 

Having shown in Figures (0) and (|5]a) the general trends of the average 
number of partons in the shower and of the average number of clans as functions 
of rapidity and virtuality variables predicted by the analytical calculations 
according to Section 4 and our previous work||], we are now ready to calculate 
the average number of partons per clan in the same variables, i.e., n c (Ay, W), 
simply by forming the ratio of the above two quantities. 

n c (Ay, W) values in rapidity intervals for W = 50 GeV (solid line), W = 
100 GeV (dashed line) and W = 500 GeV (dotted line), with parameters A 
and a equal to 2 and 1 respectively, are plotted in Figure (|5|b). The result is 
reasonable although not fully satisfactory. An anomalous behavior is clearly 
visible for large rapidity intervals (n c (Ay, W) increases as Ay — > fps in contrast 
with what one would expect from standard clan structure analysis at parton 
level) and for small rapidity intervals (n c (Ay, W) is shown to point to constant 
values, which differ for different W values in the limit Ay — > in contrast with 
the expected W independent n c (0, W) = 1 value). In order to cure the first 
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anomaly we decided to use the definition from the CPD form, Eq. (pip , which 
is relevant in fps, and accordingly to use Eq. fl23|) for the average number 
of partons per clan. The relevance of this change is shown in Figure (|6|a), 
where X'(Ay, W) is shown to diverge as Ay — > fps as expected; because of 
this behavior, however, the procedure cannot be extended to fps. It will be 
noticed in Figure (|]b) that indeed the above mentioned anomaly is removed: 
the result becomes consistent with our scheme by assuming that only clans 
with at least one parton are allowed. 

The most natural motivation of the appearence of the second anomaly 
lies in our opinion in the approximations used in the analytical calculations 
as already shown by comparing Figures @ and (f|). The approximations 
imply deviations from the exact behavior for Ay < fps. In order to check this 
guess we propose to calculate the average number of clans and the average 
number of partons per clan with the MC. Results of the MC calculations are 
shown in Figures (0a) and (0b). Their trends fully coincide with standard 
results already observed in clan structure analysis applied to parton showers 
in quark and gluon jets. The predictive power of the model under investigation 
is therefore confirmed. In addition, for small rapidity intervals the expected 
n c (Ay = 0, W) = 1 limit is recovered. 

We have shown that the parton shower model, which we built by assuming 
QCD-inspired dependence of the splitting functions in virtuality and in rapidity 
and Sudakov form factors for their normalization ("essentials of QCD"), and 
by allowing at each step in the cascading local violations of energy-momentum 
conservation law, but requiring its global validity, has an extraordinary pre- 
dictive power in regions not accessible to full perturbative QCD. In particular, 
we performed analytical calculations of the rapidity and virtuality dependence 
of the average number of clans, of the average parton multiplicity and of the 
average number of partons inside a clan in a parton shower initiated by an 
ancestor parton of given virtuality. Clans here should be intended of course in 
a more general sense than that usually referred to in the literature and linked 
hystorically to NB regularity: they are the natural language for describing a 
two step process in terms of a CPD and can be identified in our parton shower 
model with effective intermediate independent parton (gluon) sources in the 
shower. Results of our analytical calculations are consistent with what we 
know on clan properties in single quark and gluon jets disentangled by using 
a jet finding algorithm and analyzed at parton level by assuming the validity 
of GLPHD. In addition it should be pointed out that, even if the change of 
A and a parameters do not influence the qualitative behavior of the average 
number of clans and the average number of partons per clan in rapidity and 
virtuality, it actually affects their numerical values. Since A and a parameters 
control the length of the parton shower and the development of parton cascad- 
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ing inside a clan, would one describe a two parton species shower initiated by 
an ancestor gluon or quark, this flexibility of the model could be very useful 
in order to describe jets of different origin. This mechanism could be also very 
helpful in a more general framework for explaining the observed different be- 
havior of clans in different classes of high energy collisions. In summary, these 
findings suggest that the dynamics of multiparticle production is controlled 
by a QCD inspired two step process dominated by independent intermediate 
parton (gluon) sources, which we call clans, and by their cascading. The fea- 
tures which make the model appealing are its simple structure, the possibility 
to calculate in its framework analytically important physical quantities and its 
flexibility. 
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Figure 1: Integration domain in the (£,//) plane for n c (Q > 2, a) (see text, 
Eq. (|88D). Dotted lines indicate phase space domains where the integrand 
assumes different values. 
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Figure 2: Average number of partons in the shower, n(y c , W), as a function of 
the width of the rapidity interval y c obtained analytically in the GSPS model 
with A = 2, a = 1 at different maximum allowed virtualities W = 50 GeV 
(solid line), 100 GeV (dashed line) and 500 GeV (dotted line) (see Eq. (f2|)). 
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Figure 3: Same as in Figure @, but for the Monte Carlo version of the model 
at different maximum allowed virtualities W = 50 GeV (diamonds), 100 GeV 
(squares) and 500 GeV (circles). 
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Figure 4: Normalized average number of partons in the shower, 
n(y c , W)/n(ips, W) as a function of rescaled rapidity interval, y c /fps for A 
= 2, a = 1 at different maximum allowed virtualities W = 50 GeV (solid line), 
100 GeV (dashed line) and 500 GeV (dotted line). 
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Figure 5: a): average number of clans, N"(y c ,W), defined by binomial con- 
volution as in Eq. Q2"%D, as a function of the width of the rapidity interval y c 
at different maximum allowed virtualities W = 50 GeV (solid line), 100 GeV 
(dashed line) and 500 GeV (dotted line). Analytical solution with A = 2. The 
solution is of course independent of a. b): corresponding average number of 
partons per clan, n c {y c , W), as a function of the width of the rapidity interval 
y c at different maximum allowed virtualities W = 50 GeV (solid line), 100 
GeV (dashed line) and 500 GeV (dotted line). Analytical solution with A = 2 
and a = 1. 
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Figure 6: a): average number of clans, X'(y c , W), defined by requiring that the 
total MD be a CPD as in Eq. Q2"T|), as a function of the width of the rapidity 
interval y c at different maximum allowed virtualities W = 50 GeV (solid line), 
100 GeV (dashed line) and 500 GeV (dotted line). Analytical solution with A 
= 2. The solution is of course independent of a. b): corresponding average 
number of partons per clan, n c (y c , W), as a function of the width of the rapidity 
interval y c at different maximum allowed virtualities W = 50 GeV (solid line), 
100 GeV (dashed line) and 500 GeV (dotted line). Analytical solution with A 
= 2 and a — 1. 
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Figure 7: Same as in Figure (|6]) , but Monte Carlo results at different maximum 
allowed virtualities W = 50 GeV (diamonds), 100 GeV (squares) and 500 GeV 
(circles). 
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